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Abstract. The shift operator for a quantum lattice current algebra 
associated with si (2) is produced in the form of product of local fac- 



■ tors. This gives a natural deformation of the Sugawara construction 
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, for discrete space-time. 

^ ■ Introduction 



The Current Algebra provides the chiral dynamical variables for a generic con- 
formal field theory model called WZNW model. Its lattice analogue, due to 
Semenov-Tian-Shansky, proved useful for elucidating the quantum group struc- 
ture in this model [AFSV]. In the subsequent papers [AFS, FG, BC, AFFS] 
some general properties of this algebra and its representations were discussed. 
However, these considerations covered kinematical aspects of the lattice model 
while such a basic dynamical object as the hamiltonian density remained un- 
available. Here we address this problem making use of our experience in a 
simpler abelian case [FV93, V95]. Following the general philosophy worked 
out in these papers we construct a spatial translation operator W which simul- 
taneously generates the temporal shift. We find W to be a product of local 
factors over the lattice. This may be regarded as a multiplicative analogue of 
the Sugawara construction. 

For simplicity we confine ourselves to the simplest case of the sl{2) algebra. 
In Section 1 we recall the basic facts about the current algebra in its classical 
continuous form. Then we embed the Sugawara hamiltonian into the hierar- 
chies of conservation laws of two major integrable models which are mKdV and 
NLS equations [FT]. To make a smoother transition to the quantum case we 
present in Section 3 the classical lattice deformation of the current algebra. In 
particular, we produce relevant integrable hierarchies. The quantum case is 
treated in Section 4. 
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1 Classical Model 

The generators j a (x) of current algebra are associated with a given simple Lie 
algebra g with index a labeling the linear basis in g and the variable x running 
through the unit circle. Let and K ab be the structure constants and the 
Killing tensor of g. The defining Poisson bracket is 

if(x)J b (y)} = lftf{x)5{x -y) + 7 K ab 5\x - y) . 

The real 'coupling constant' 7 is irrelevant in classical case but comes into play 
under quantization. 
The hamiltonian 

H = ±- K ab f(x)j\x)dx 
*H Jo 

leads to a free equation of motion 

d t f(x) = {H,f(x)} = d x f(x) 

which reflects the conformal invariance in hyperbolic language. The hamiltonian 
density 

T(x) = K ab j a {x)j\x) 

is quadratic in the generators and is often referred to as the Sugawara construc- 
tion. In this paper we shall consider g to be a real form sl(2, R) of the algebra 
sl(2). 

Thus, a takes the values 3,+,— and all the functions j a are real. The 
Poisson bracket is given by 

{f(x),j\y)} = 1 5'{x-y) 

{f{x),j ± (y)} = ± 1 j ± 5(x-y) 
{j + (x),r{y)} = 2j(j 3 5(x -y) + S'(x - y)) 

{j ± (x),j ± (y)} = o 

and 

It is also useful to combine the currents into a 2 by 2 matrix 

and write the Poisson bracket in the form 

{} (x), J (y)} = \[C,J Or)- J (y)]5(x - y) + 7 C«5'(x - y) 

where the standard notation of the i?-matrix formalism is employed [FT], and 
C is the Casimir element. 
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2 Separation of Variables and Yang-Baxterization 



The above bracket and hamiltonian allow for separation of variables. Indeed, 
one may put the matrix J into the form 



with a diagonal matrix f2 solving the equation 

d x n = fa 3 n. 

The Poisson bracket for the new set of dynamical variables j = j' 3 , p, q proves 
to be 

{p{x),p{y)} = -27sign(x - y)p{x)p{y) 
{q(x),q(y)} = -27sign(a; - y)q{x)q(y) 
{p(x),q(y)} = 27(sign(a; - y)p{x)p{y) + S'(x - y)) 



{j(x),j(y)} = n/6'(x - y) 
{j(x),p(y)} = {j(x),g(y)} = o. 

while the hamiltonian density becomes 

T = j 2 +pq. 

Thus, the pair p, q completely separates from j. The p-q bracket is known to 
belong to the hierarchy of Poisson structures associated with the NLS equation 
while the density pq is a member (the momentum density) of the corresponding 
family of densities of local conservation laws [FT]. On the other hand, the j- 
bracket and the density j 2 come from the hierarchy of the mKdV equation. 
Thus, we see where sl(2) current algebra and Sugawara hamiltonian fit into the 
general pattern of Soliton Theory: 

"^WZNW = ^mKdV + P NLS' 

This will prove useful for our approach to quantization. 

The lattice formalism for the mKdV part, which is nothing but the abelian 
current algebra, was already developed in [FV93, V92, V95]. In this paper we 
perform a similar treatment of the NLS part. 

In Soliton Theory the densities of conservation laws come from the asymp- 
totic expansion of the trace of the monodromy matrix of the auxiliary linear 
problem. For the NLS equation this problem reads 
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The matrix f2 being diagonal, this auxiliary problem is gauge equivalent to 

(d x + J + Xa 3 ) * = 0. 

Thus, we see that NLS part -f^g °f the Sugawara hamiltonian is provided by 
'Yang-Baxterization' of the current 

J ~> J + A03. 

In the next Section we shall do the same on the lattice. 

3 Lattice Model 

We discretize the circle introducing the spacial variable taking integer values 
running from 1 to N. The real dynamical variables will be denoted by a n+ i , (3 n 
with integer n; it is understood that 

a n+N+ i = a n+ i 

One may say that integers label vertices while half-integers stand for edges. Or 
vice versa. One reason for using different notations for dynamical variables with 
integer and half-integer subscripts is mere convenience which becomes evident 
when the Poisson bracket is displayed: 

K-±> a n+i} = -27a n _i«„ + i 
{a n _i,0 n } = -2 7 a n _i/3 n 
{I3 n ,a n+ i} = -2 7 /3„a n+ i 

All brackets not listed are zero. It is clear that every variable has nontrivial 
brackets only with the two neighbours in either direction. 

One can recognise here the so called Flaschka variables for the Toda model. 
However, the hierarchy we will deal with is different from that of the Toda 
equations. 

To see what this lot has to do with the Current Algebra we arrange dynam- 
ical variables in two matrices 




( 1 M 
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The Poisson relations for them 



12 12 
{B2n, B2n} = 7[ r 12, B2nB2n] 

12 12 
{C2n-1, C2n-l} = 7[ r 21 , C2n-lC2n-l] 

12 12 
{B2n, C2n-l} = 7-B2ri?*12C2n-l 

12 12 
{C2n+l,B2n} = lC2n+\T2\B2n 

employ the major ingredient of g-deformations, namely the classical r-matrices 



ri2 











°\ 





1 

2 


2 











l 

2 





V o 








\) 



where P is a permutation. 
The product 

satisfies the Poisson brackets 



T21 = P\2r\2P\2, 



Jn — B2 n C2n-l 



12 12 12 

{Jn, Jn} = l{r\2JnJ2n ~ JnJ2n?"2l) 

12 12 
{Jn+1, Jn} = lJn+lT2lJn 

which turn into the Current Algebra in the continuum limit 

J n ~ I + AJ(x). 

This is what usually is called the Lattice Current Algebra. However, it is not 
clear whether one gains anything reducing B-C-algebra to the J-one. This time 
we prefer to deal with somewhat more transparent -B-C-algebra but we could 
do with the J-one instead. 

To produce relevant conservation laws we introduce the so called 'transfer- 
matrix' 

t{u) =trnr 3 52n*7~ a3 C 2 „-i 

n 

with 'spectral parameter' to entering in £ and rj in such a way that 

e + rj 2 = 2 



i 

V 



U). 



It turns out that 

(i) t(co) is a Poisson commuting family: 

{t(cj),t(co')} = 0, 
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(ii) in the continuous limit it turns into the trace of monodromy matrix of 
the continuous auxiliary linear problem of Section 1 provided 



to ~ 1 + AA 



(iii) it is a power series in u 



N/2 
-N/2 



2k 



with 




(ii) is obvious, (iii) is almost so, (i) can be verified along the guidelines of [FM]. 
We shall not go into further details because the model in question actually 
belongs to the same hierarchy as the Ablowitz-Ladik's model [SV]. 
In the continuous limit we have 



as should be expected. 

We have obtained the hamiltonian of the classical lattice model which plays 
the role of the NLS part of the Sugawara construction for the lattice current 
algebra. The corresponding mKdV part can be found in [V92]. However, the 
equations of motion produced by these hamiltonians are quite complicated and 
turn into simple free equations only in the continuous limit. It was realized in 
[FV93] that the discrete time equation 



is a better option. In other words, the discretizing of space should be accom- 
panied by the discretizing of time. The last equation is especially transparent 
in the quantum theory where the spacial shift operator W such that 



We shall find this operator in the next Section. The expression for the classical 
lattice hamiltonian will prove to be a useful hint in our search. 




NLS 



Jn{t + A) = J n+ i(t) 



W' 1 J n W = J n+1 



is taken to define the time shift as well 



J n (t + A) = W- 1 J n {t)W 



Jn(0) = J, 
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4 Shift Operator 

The quantum lattice current algebra inherits the notation a-(3 for generators 
together with the way they are enumerated while the Poisson relations turn 
into their most natural quantum counterparts 

[0n-nA»] = 0?-g - V n _I , 

with the deformation parameter q combining the coupling constant 7 and the 
Planck constant H in the usual way 

q = e in ~<. 

The consistency of these commutation relations becomes more apparent as soon 
as one rewrites them in ii-matrix form 

12 21 

R\2B2nB2n = B2nB2nR\2 

12 2 1 

-R2lC2n-lC2n-l = C2n-lC2n-l-R21 

2 11 2 

C2n-\B2n = B2nRl2C2n-l 

2 1 1 2 

-B2nC2n+l = C2n+\R2\B2n 

where matrices B,C are built of a,/3's in literally the same way as in the classical 
case of Section 2. The i?-matrix involved is the sl(2) one 



R12 



( 



V 



_ 1 
q 2 



I _3 

p — q 2 
_ 1 
q 2 



q2 j 



and it is needless to say that the associativity of the B-C-algebra is due to the 
Yang-Baxter equation 

R12R13R23 = R23R13R12 

fulfilled by R. 

The way variables separate in the continuous model and the belief that 
integrals of local densities on the lattice turn into products of local factors 
suggest that the shift operator 

a n _iW = Wa n+3 



(3 n _ lW = Wf3 n+1 
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should decompose into a product of two commuting factors 



W = UV = VU 

of the form 
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It only remains to find suitable functions and a. It is shown in Appendix B 
that the functions doing the job are 

M , (log(-^)) 2 
6(z) = exp ■ 



41ogg 
and 



1 f a 
a(z) = exp - / 

*± J — c 



sinh(vrC) sinh(77iC) C 
where the contour of integration rounds the singularity at £ = from above. 

Let us conclude with some remarks. It must be said that although the 
solution looks like a straightforward remake of the one for the £7(1) case [FV93] 
it is not. The logics of the abelian case would rather favour another local 
decomposition for the shift operator, with functions 9(— z) and cr(— z) instead 
of 9{z) and <r(z). Of course, the minus in the argument of a is a lot more 
important than it may look. In particular, the r.h.s. of the functional equation 
fulfilled by a(z) (see Appendix) comes directly from the density of the classical 
lattice hamiltonian 



1 + 1 2 + 0(3/a 

This correspondence principle plays a major role in the detailed study of the 
classical limit which will be presented elsewhere. 

As we said, one could do with the J-picture from the very beginning. This 
would eventually lead to the following decomposition 

for the J-shifting operator 

J n W = WJn+i. 

It is easy to check that all entries of this decomposition do belong to the J- 
subalgebra of the i?-C-algebra. However, it is not that apparent why the re- 
duction B-C — ► J should eliminate some factors in W and modify some of 
remaining ones. Subtle things like central elements are partly responsible for 
this. Which in turn has a lot to do with another important issue which we 
deliberately left aside. Indeed, we did not pay proper attention to the peri- 
odic boundary condition. What was then meant by ordered products over the 
lattice? For the answer to this question we refer the reader to [V95]. 
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Appendix. Full Pull 

In this paper, as well as on a few earlier occasions [FV93, FV94, F, FV95], we 
relied on certain '(/-special-functions' satisfying functional equations of the form 

f{qz) -d{z). 



In particular, the functions 9 and a of the last Section were given by d(z) = z 
and d(z) = j-^ respectively. Although the issues of the solvability and the good 
choice of solution are important in their own right, for practical purposes (read 
formal computations) one seldom needs anything else than just the equation 
itself. Let us show how it goes in our case. 

Due to the locality of commutation relations the equations 



a n _iW = Wa n+ 3 
P n _ lW = W(3 n+1 



easily reduce to 



and 

The mere form of the functional equation suggests that any computation ought 
to be a sequence of just two elementary steps: 

(i) wf(z) = f(z) w if wz = zw, 

(ii) xf(y) = f(y) d{qy)x if x, y make a Weyl pair: xy = q 2 yx. 

Since any two of the a's either commute or make a Weyl pair, the first trans- 
lation comes easy: 

a n _ h e{qa n+3 _a-l h )e{qa n+1 _a-\) 

The second one is more tricky. We cannot pull (3 n _ 1 through a(t n+ i) 
straight away because /3 n _ 1 and t n+ i neither commute nor make a Weyl pair. 
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Nevertheless, we have a good supply of operators making 'good' pairs with both 
t„_i and i_ , 1 which, by the way, between themselves are a Q-oscillatorf] 

III (-) fit I r} I I 



qt n+ ^t n _,-q \_it n+ i =q-q 

Among them are: 

(i) all the a's 

[a,t] =0, 

(ii) the /? which is 'between' them 

P n t n+ i = q\ + i/3 n , 

(iii) another operator c n = q(t n+ it n _i — 1) 

cj n+ i = q\ + ic n 
which is a familiar satellite of (/-oscillators. 
So, we express j3 n _ l via 'good' operators 

and get 

/3 n _iO-(t n+ i) 

2y~2 2 2 2 2 

= a(t n+ i ) {{(i~X-i ~ t n+ i_p- l )a n _ h )) . 

Similarly, 

^(t n _i)f3 n+1 = (a n+ i(-p-%_i +t nH (3~ 1 )) a(t n _i). 

The rest 

is a variation on the same theme. 

Unfortunately, all this looks like a hat-trick rather than conscious approach. 
A more systematic paper making better use of YangBaxterization is on the 
authors' agenda. 

"The remaining nontrivial commutation relations governing the algebra of t's, those for 
the neighbours twice removed, seldom participate in computations. Their role may be seen 
in taking care of the associativity of the algebra. Anyway, their explicit form can be found in 
[V92]. This algebra is sometimes referred to as the Lattice Virasoro Algebra for in a certain 
continuous limit, different from the one of the present paper, it turns into the Virasoro algebra 
with a nonzero central charge. 
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